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THE FATE OF C0 2 FROM THE NATUNA GAS PROJECT IF DISPOSED BY SUBSEA SPARGING 


We investigate the consequences of disposing C0 2 waste gas from the Natuna project by 
sparging C0 2 into the ocean at 140 m depth near the site. The large reservoir of natural gas 
discovered near Natuna Island in the South China Sea contains over 70 percent C0 2 . Proposed 
levels of production amount to about 0.4 percent of the current total global C0 2 emissions. This 
would make Natuna the world's largest point source emitter of C0 2 and raises concern for the 
possible incremental impact of Natuna on the C0 2 greenhouse problem. While the base case produc- 
tion scheme calls for atmospheric venting of C0 2 , an alternative has been proposed in which C0 2 
would be sparged into seawater so that, perhaps, C0 2 would never appear in the atmosphere. 

Here we develop models for the chemistry and transport of C0 2 in seawater. The 
chemistry is well understood and has been investigated extensively in connection with the C0 2 
Greenhouse issue. Transport models are more complex because currents can carry C0 2 over the 
entire basin of the South China Sea within a period of months. Data for the three dimensional, 
time dependent current distribution do not exist. Our models consider diffusion and advection 
with parameteri zat i ons characteristic of expected dynamics. While simplified, the models are 
readily interpretable, and serve to estimate the order of magnitude of likely transport effects. 

Carbon in seawater is not inert, C0 2 rapidly partitions itself into carbonate and 
bicarbonate species, changing the oceanic pH. If the pH falls below 7.4, calcite in shells 
begins to dissolve. Increases in carbon become magnified in their effect on the partial pressure 
of dissolved C0 2 . When the partial pressure of C0 2 in the ocean exceeds the partial pressure of 
C0 2 in the air, C0 2 degasses to the atmosphere. As sparged C0 2 builds up in the basin it will 
begin to degas. After a time which we refer to as the "retention time" C0 2 will degas at a rate 
equal to the injection rate from the sparger. 

The principal conclusions from our models are that (1) the retention time is only about 
ten years or less, and (2) addition of C0 2 raises the acidity of seawater sufficiently to cause 
dissolution of calcite over an area of order 1000 square km in size. These effects indicate that 
sparging offers no advantage over direct atmospheric venting of C0 2 . 

At this time, no further work on this project is planned. 


Brian P. Flannery, Andrew J. Callegari, Bahlin Nair, and Wayne G. Roberge 
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1. INTRODUCTION 

Initial plans for development of natural gas discovered off Natuna 
Island in the South China Sea produce as waste gas approximately 65 x 10 6 
metric tons of CO 2 per year (1.5 x 10^ moles CO 2 per year). Release of this 
amount of CO 2 to the atmosphere raises concern with respect to its effect on 
the CO 2 greenhouse problem. Global fossil fuel emissions of CO 2 currently 
amount to about 1.8 x 10^ metric tons per year, so Natuna would represent 
about 0.4% of world CO 2 emissions, but would occur as a localized point 
source. 

Various descriptions of the Natuna project cite different values for 
the amount of CO 2 produced. Slight differences arise depending on the assumed 
production capabilities of the offshore platforms, major differences arise 
depending on the number of platforms in use. Our figure assumes total 
production of 5540 million standard cubic feet per day of raw gas from three 
platforms required to handle the "Pipeline" and "LNG" projects. None of the 
major conclusions of this report are materially affected by slight changes in 
the overall production rate. 

A possible disposal option that might not lead to atmospheric 
release of CO 2 involves pumping waste gas back to the sea floor. The idea of 
using the deep ocean to dispose of industrial CO 2 was proposed by Marchetti 
(1977). Hoffert et al . (1979) investigated the proposal in more detail and 
found that injection of CO 2 onto the seafloor at 4000 m depth would prevent 
CO 2 from entering the atmosphere for a time of order 1000 years. In this 
study, we investigate the consequences of CO 2 addition to seawater in the 
shallow basin near the Natuna site. We focus on two principal features (1) 
the role CO 2 addition on seawater chemistry, especially pH, and (2) the 
retention time for CO 2 to remain dissolved in seawater before degassing to the 


atmosphere. 
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Figure 1 shows geographical features of the site in the South China 
Sea basin. Over much of the basin the seafloor is not at great depth, but 
lies on a flat, shallow shelf where the water depth does not exceed 200 m. 
Monsoons sweep the basin at six month intervals driving and seasonally 
reversing surface currents that circulate in a gyre over the basin (Wyrtki, 
1961). We consider discharging CO 2 rich waste gas into seawater near the 
seafloor through a sparger system as illustrated in Figure 2. The sparger 
itself consists of a ring approximately 10 km in diameter from which 1.5 x 
10^ moles per year of CO 2 would be degassed at a depth of 140 m. Vertically, 
there are three regions of interest. The sparged CO 2 rises approximately 30 m 
in a bubble plume before dissolving. At the surface wind and wave action 
create a well mixed layer of about 30 m depth. Finally, in the middle regions 
there exists a stably stratified region of about 80 m thickness. Currents 
near the sparger vary in magnitude and direction with height and time of 
year. Typical values at depth are of order 10-50 cm s“*. At the bottom 
currents are persistently headed southwest, toward Borneo 500 km distant. At 
the characteristic velocity of 30 cm s"* flow traverses the basin in only 
20 days. Thus, advective transport can distribute material injected at Natuna 
over a basin scale in a time of order months. 

Adding CO 2 to seawater decreases the pH and increases the partial 
pressure P(C02) of dissolved CO 2 . As described in Section 2, the chemistry of 
CO 2 in seawater is well understood, and has received much attention recently 
in the context of the role of the oceans as a sink for atmospheric CO 2 in the 
CO 2 greenhouse problem. At the ocean's surface, gaseous CO 2 is exchanged with 
the atmosphere whenever a difference exists between the partial pressure of 
C0 2 in the atmosphere P A (CO 2 ) and the partial pressure of CO 2 dissolved in 
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surface seawater P$ (CO2). CO2 flows from the reservoir with higher pressure 
to the reservoir with lower pressure. Thus, CO2 will degas to the atmosphere 
when P s (C0 2 ) >P/^ (C0 2 ) • 

In Section 3 we define a CO2 retention time for the Natuna problem 
according to the concept of time taken for seawater to reach a steady-state 
concentration of CO2. For Natuna that occurs when the ocean has taken up so 
much CO2 that it degasses to the atmosphere at the same rate that spargers 
inject CO2 at the bottom. Estimates in Section 3 illustrate the central idea 
of retention time, and establish its order of magnitude for a seawater basin 
in which CO2 rapidly mixes throughout the entire volume. 

A complete investigation of mixing and transport would require 
detailed current information in three dimensions and time throughout the South 
China Sea basin. Such data does not exist, and, even if it did, to use it 
would require a massive computational simulation. Instead we develop insight 
and understanding using simplified models and analysis. Nonetheless, these 
simpler models allow us to estimate the magnitude of the retention time and 
impact on pH. In Section 4 we consider several classes of models that include 
the effects of transport by diffusion and advection. Typically, such models 
predict retention times longer than estimates based on the simple well -mixed 
model described in 3, but retention time still remains short compared with the 
timescale of climate change in the CO2 greenhouse problem. On the other hand, 
transport models show greater pH reductions near the sparger than the well- 
mixed model described in Section 3. 

Finally, in Section 5 we discuss these results in the context of the 
global CO2 greenhouse problem. Our conclusion is that atmospheric discharge 
is preferable to seawater sparging. In any event, seawater sparging is nearly 
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equivalent to atmospheric release since the retention time is short compared 
with times relevant to the buildup of atmospheric C0 2 influencing climate 
through a greenhouse effect. 

2. C0 2 CHEMISTRY IN SEAWATER 

The chemistry of C0 2 in seawater has been discussed at length in 
many recent references. Our discussion follows that of Takahashi et al . 

(1980), also see Baes (1982) who developed a graphical representation of the 
results that is particularly convenient for us. For the Natuna project, the 
salient features of C0 2 chemistry are the variation in pH and P(C0 2 ) produced 
by changes in the total amount of carbon, TC, dissolved in seawater. Changes 
in pH directly influence chemical reactions, such as the oxidation of H 2 S and 
solubility of calcium carbonate, that influence toxicity. P(C0 2 ) in surface 
water controls the direction and rate of gaseous C0 2 exchange between air and 
sea that determine the length of time that sparged C0 2 will remain in the 
ocean. The essential result from consideration of C0 2 chemistry is that small 
increases in TC cause decreases in pH and large increases in P(C0 2 ). 

C0 2 added to seawater rapidly dissolves to form hydrated C0 2 and 
carbonic acid H 2 CO 3 . In turn, the carbonic acid can dissociate to form bicar- 


bonate HCO3 and carbonate CO^" ions 


C0 2 (g) + H 2 0 l H 2 C0 3 (aq) 

( 2 . 1 ) 

H 2 C0 3 (aq) * H + + HCO 3 

( 2 . 2 ) 

hco: * H + + C0?‘ 

3 «- 3 

(2.3) 
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where (g) and (aq) denote gaseous and dissolved aqueous CO 2 . All reactions 
proceed rapidly enough that the ultimate distribution of carbon among the 
three species can be described by equilibrium constants 


K x = [H 2 C0 3 ]/P(C0 2 ) (2.4) 

K 2 = H + [HC05]/[H 2 C0 3 ] (2.5) 

K 3 = [H + ][C0§-]/[HC0 3 ] (2.6) 

Note that reactions (2.2) and (2.3) depend directly on the oceanic pH. The 
equilibrium "constants" depend on temperature, pressure, and salinity. 

For a given addition of C0 2 to seawater we want to determine the 
resultant change in pH and the partial pressure P(C0 2 ). In general, we can 
predict those quantities given the titration alkalinity TA, the temperature, 
and the total amount of carbon TC, defined as 

TC = [H 2 C0 3 ] + [HCO^] + [C0§’] (2.7) 

TC measures carbon concentration per unit mass, occasionally we also reference 
carbon abundance with respect to volume, tc = pTC, where the density of sea- 
water is 1.025 gm cm - ^. At Natuna a typical value for TC is 1930 pmole kg'* 
or 1.98 moles m - ^. Titration alkalinity, a quantity measurable with high 
precision, represents the net molar concentration of positive ions whose 
abundance is not sensitive to pH. The excess cationic charge resulting from 
dissociation of strong electrolytes in seawater is balanced by the anionic 
charges which are mainly generated by dissociation of carbonic and boric 
acid. Thus, the titration alkalinity is: 
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TA = [HCO 5 ] + 2 [C 0 §'] + [H 2 B 05 ] + [OH - ] - [H + ] ( 2 . 8 ) 

Effects of other minor species such as silicic and phosphoric acids are 
neglected. In seawater TA is dominated by dissociated salt species. Conse- 
quently, TA is a function of salinity, and does not depend on TC. Given TC, 
TA, temperature, and functional relations defining the equilibrium coeffici- 
ents, the set of equations can be solved for consistent values of pH, P(C0 2 ), 
and the abundances of the carbon species. For the equilibrium coefficients we 
use recent values supplied by Taro Takahashi of Lamont Doherty (Takahashi 
et al., 1982). 

Results for effects of local chemistry are shown in Figure 3-4 which 
illustrate the dependence of pH and P(C 02 ) on TA and TC for temperatures of 

20, and 30°C, representative of surface and deep water near the Natuna 

drilling site. Based on the global correlation between salinity and TA 
(Takahashi et al., 1982), Takahashi estimates that TA near Natuna should be 
about 2300 yeq/kg. Air sea exchange tends to maintain equatorial waters 
nearly in equilibrium with atmospheric P(C 02 ) at a value of 340 yatm. Thus, 
surface waters at Natuna have a value TC = 1930 y mol /kg as shown in Figure 3. 

The important feature shown in Figures 3-4 is the strong increase in 

P(C 02 ) and decrease in pH for increases in TC at fixed salinity and tempera- 
ture which arises from the well known buffer factor in seawater. Notice that 
increases of only 10% in TC increase P(C 02 ) from 340 yatm to over 1000 yatm 
and decrease pH from 8.2 to 7.4. In the stipled regions of Figures 3-4 
calcium carbonate is undersaturated so that shell material will dissolve, 
unless biological activity offsets the dissolution rate. The figures illus- 
trate that CO 2 is less soluble at higher temperature. Thus, for a given 
absolute value TC, pH is lower and P(C0 2 ) higher at higher temperature. 
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That P(C0 2 ) increases rapidly with TC is well known (e.g. Weiss, 


1974); the ratio of proportional variation for small increments is known as 
the Revel le factor y, defined as 


AP(C0 2 )/P(C0 2 ) = y ATC/TC 


(2.9) 


For surface waters at Natuna we calculate y = 9.2. For surface water, the 
Revel le factor ranges from 8 in warm waters to 15 in cold waters, globally 
averaging 10 (Broecker et al., 1979; Takahashi et al., 1980). The Revelle 
factor is especially useful for analysis of the C0 2 retention time as 
described in the next section. 


We can use results of this section to estimate the impact of 


sparging in the near vicinity of the sparger itself, on length scales of order 
10 km when local seawater transport is dominated by advection. We assume that 


the C0 2 influx, Q = 1.5 x 10^ moles/year, is released uniformly into a sea- 


water flowing over the sparger 10 km wide by 30 m high, the height of the 
bubble plume (hg). Then the change in tc caused by sparging is given by 


CO^ flow rate Q 


( 2 . 10 ) 


water flow rate D*h D .v 


-1 


.3 

Ate (moles m ) = 1.6 ( 



-1 


Figure 5 shows tc vs. velocity in the range 1-100 cm s , and the correspond- 


ing variation in local P(C0 2 ) and pH. Note that if local stagnation events 
occur they will be accompanied by large changes in local pH. Onsite measure- 


ments indicate that velocities less than 10 cm s"* do occur occasionally. 


t 
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3. THE SEA-AIR EXCHANGE OF C0 2 AND THE C0 2 RETENTION TIME 

Gaseous C0 2 migrates across the sea air surface interface at a net 
rate F SA proportional to the difference between the partial pressures of C0 2 
in the two reservoirs 



p s (C0 2 ) - Pfl(C0 2 ) 

340 yatm 


■] 


(3.1) 


The expression has been scaled to units appropriate for the annual mean yearly 
average partial pressure of atmospheric C0 2 , P A (C0 2 ) = 340 yatm. Empirical 
calibrations give E = 20 moles/(m 2 yr) (Peng et al . 1979). For surface condi- 
tions at Natuna, tc eq « 2.0 moles m" 3 is required for P s (C0 2 ) = 340 yatm, at 
which F^ a vanishes. Surface waters over most of the worlds oceans maintain a 
carbon abundance such that F$ A « 0. 

To estimate the appropriate timescale, t eq , for air-sea gas exchange 
to equilibrate, consider a situation in which a well-mixed layer of seawater 
of depth h contains no C0 2 initially, and is in contact with air at P A (C0 2 ) = 
340 yatm. For typical oceanic conditions the well-mixed layer has a depth h = 
75 meters. To reach equilibrium the mixed layer requires the addition of 
h*tc = 140 moles m" 2 , and the initial filling rate is 20 moles/(m 2 yr). Thus, 
a characteristic estimate for the equilibration time scale is 


h tc 


eq 


eq 


SA 


7 yr 


(3.2) 


After time of order t pa the net exchange of C0 2 between air and sea vanishes 

CV-j 

and tc reaches steady state in the ocean. 
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While the estimate for air sea exchange equilibration time in (3.2) 
considered no additional sources, for the Natuna problem the CO 2 concentration 
in seawater will not reach equilibrium until the ocean waters lose CO 2 at a 
rate identical to the influx from the sparger Q. In later sections we will 
consider more detailed models including diffusion and transport, our goal here 
is to obtain simple estimates for the characteristic timescale of the problem. 
We estimate the retention time tr for CO 2 by first estimating what the carbon 
concentration tc e q in surface waters must be in order for surface outflux to 
equal Q. We then estimate the equilibration time assuming by that the neces- 
sary increase in tc occurs at a rate Q. 

Nothing in the estimates used here actually allows us to determine 
the basin size, or surface area over which degassing occurs, we will attempt 
to estimate basin size in later sections. Here we solve for retention time as 
a function of assumed surface area for a circular region of arbitrary radius 
R. Because the carbon chemistry of seawater produces such a strong increase 
in P(C 02 ) with increase in total CO 2 the retention time is very short for 
small basin sizes. For large basin size, the retention time approaches a 
limiting value that does not depend on radius. Thus, we will obtain useful 
(lower) bounds on the retention time even without knowing the size of the 
affected region. 

For a circular surface area of radius R the net outflux of CO 2 from 
sea to air will equal the sparger influx for 


ap(co 2 ) 

340 yatm 


„ R f sa - q 


2.4 


l 1.5 X 10 12 moles yr -1 


100 km 
R 


(3.3) 


10 - 


? 


To estimate the increase in TC required to produce the change aP(C 02) we must 
solve the non-linear chemistry equations of the previous section. However, 
for small changes in TC we can use the Revel le factor to determine the 
concentration change directly. 


n 


R 2 E 


(tt) ■« 


for 


ATC 

T (T 


« 1 


AtC 

tc 


0.26 


1.5 x 10 12 mole yr" 1 


100 km 


(3.4) 


Note that initially P S (C02) = P/\(C02) = 340 uatm, and that Atc/tc = aTC/TC. 

Given the increase in TC, either from (3.4) or by solution of the 
nonlinear equation (3.3) we can crudely estimate the time for the ocean to 
reach equilibrium as in the discussion for t e q in (3.2). Assume that the 
entire basin of volume nR^h must increase its CO 2 concentration by Ate, and 
use the initial influx rate Q to estimate the filling time 


Tr> = 


uR h Ate 


(3.5) 


For the Natuna basin we use h = 140 m. If the required change Ate is small, 
then we can use (3.4), based on the Revelle factor, to determine tr. 


h tc 
t R y e 


1.4 yr 


for 


Ate 

tc 


« 1 


(3.6) 


Notice that this estimate is independent of the basin size R and the influx 
rate Q. 
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The crude estimate derived here represents an approximation corres- 
ponding to adding CO 2 to the entire cylindrical volume of the basin as though 
the concentration remained well mixed at all times. This is unlikely to occur 
both because the water is stably stratified vertically and because advective 
transport controls horizontal distribution. In later sections we consider 
such additional effects. Nonetheless, these estimates place the character- 
istic timescale in perspective, it is of order years. 

Using these estimates, and accurate solution of the nonlinear equa- 
tion (3.3) to determine tc appropriate to equilibrium degassing for P(C02), we 
illustrate several features of the results in Figure 6. In Figure 6a and b we 
show the required change in P(C 02 ) and tc as a function of basin size. 

Figure 6c shows the retention time estimated from (3.5). Given the increase 
in tc we also show the pH corresponding to equilibrium conditions. The 
important point is that, if rapid mixing occurs on small length scales, of 
order of sparger diameter, 10 km, then the degassing time is very short 
regardless of the exact size of the affected region. Over larger scales the 
degassing time will still be only a few years if CO 2 reaches the surface. 

Detailed models of transport phenomenon as discussed below increase 
these time estimates. 

4. TRANSPORT MODELS 

Injected CO 2 from the sparger dissolves into seawater in a bubble 
plume within 30 m of the seafloor. Measurements of vertical profiles of 
temperature and density at the site indicate that below a well-mixed surface 
layer 30 m deep the water is stably stratified, so vertical mixing must 
proceed by eddy diffusion, sporadic upwelling events, or complex non-local 
transport by currents. Our estimates below will show that vertical diffusion 
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times of order a few years are required to bring CO 2 to the surface. As 
described in the introduction, measured horizontal velocities of order 
30 cm s _1 carry material over basin scale distances in only months. Site 
measurements also indicate that current varies significantly with depth and 
time. Consequently, a realistic model of seawater transport of C0 2 in the 
South China Sea requires current information covering the entire basin in 
three spatial dimensions and time. 

To date no complete data set exists to carry out a realistic trans- 
port calculation on the scale of the basin, and such a calculation presents a 
formidable numerical exercise even if the data were available. To gain 
insight into the likely retention time for C0 2 we investigate a hierarchy of 
simpler models based on assumed current and diffusion profiles, that do 
include time dependence. Cumulatively, they indicate that C0 2 retention times 
will be of order a decade at most, barring unlikely current patterns, or 
anomalously low vertical diffusion rates, as we describe. 

The models described below have analytical solutions for the time 
dependent variation of concentration when can be approximated by a linear 
dependence on concentration, using the Revelle factor in (3.1). However, 
numerical solutions are simple to obtain, even with non-linear relations for 
the loss term, and simple to interpret based on known characteristic behavior 
for the diffusion equation. Table 1 contains a listing of model parameters 
and typical values for their magnitude. 

K( V) and K(H) representing vertical and horizontal diffusion in the 
ocean are the only parameters listed in Table 1 that have not yet been dis- 
cussed. Below the mixed layer of the ocean, variation of temperature and 
density produce stable stratification that inhibits vertical flow of any sort, 
including diffusive mixing. K(V) is particularly important in this problem 
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since it sets the timescale for carbon to diffuse to the surface where losses 
can occur by air-sea exchange. The value we list, 4000 m 2 yr“^ 

(- 1.3 cm s ) is typical of vertical diffusion through stably stratified 
layers in the ocean as a whole (Broecker and Peng, 1982). A recent study of 
Berelson et al . (1982), based on the vertical distribution of radon-222 in the 
Santa Barbara basin of the California coast, found K(V) = 4 cm 2 s' 1 at water 
depths of 500 m. No measurements of K(V) exist for the Natuna basin. 

It is possible that wave action and currents in a confined basin 
make K(V) become larger, reducing the upward diffusion time, or that special 
circumstances cause K(V) to be smaller than for the ocean as a whole. No a 
priori theory exists that can accurately predict values for K(V). 

Horizontal diffusion is known to be dramatically larger than 
vertical, K(H) = 10^ K(V), is generally accepted (Kuo and Veronis, 1973; 
Sarmiento et al., 1982). However, the model results are less sensitive to 
K(H), since advective current flow tends to dominate diffusion as a cause for 
horizontal mixing. 


a * Models With Rapid Horizontal Mixing and Vertical Diffusion 

In the previous section we estimated the retention time of a well- 
mixed cylindrical model based on the required equilibrium concentration and 
initial influx rate. Here we consider the time history of such a model, and 
add effects due to vertical diffusion. Since surface losses rise as CO 2 
increases in the reservoir, the actual history of outgassing causes the reten- 
tion time to increase over that estimated in the previous section. 
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In this model we assume that horizontal mixing throughout the basin 
of radius R occurs rapidly at any depth, but that vertical transport occurs 
only by diffusion in the stratified layers below the surface mixed layer of 
depth h M . The C0 2 injection rate per unit volume q(z) is taken to be uniform 
per unit volume at depths within the bubble plume, z < h B = 30 m. The 
equation defining the time history of concentration is 

fr- ■ q (z) * I? K( v) !f c (4 - 1) 


where 


q(z) = — for z < h B 
irR n g 

q(z) =0 for z > h g 


Here z represents height above the seafloor in a water column of total depth 
H. Boundary conditions for (4.1) correspond to separate treatment of the C0 2 
balance in the mixed layer and zero diffusive flux across the seafloor. 


8tc(M) . 

at 


- K(V) 


ate , 

az m-h M 



(4.2a) 


ate , 
az z 



(4.2b) 


Several features of the model should be noted. First, from the form 

o 2 

of (4.1-2) the solution depends only on the basin area ttRS and volume ttR H, 
but not on the assumed cylindrical geometry; we use R simply as a convenient 
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measure of basin size. Second, in the limit where K(V) becomes very large 
this model becomes vertically (and horizontally) well mixed, as in the simpler 
model of Section 3. Mathematically, the vertically well mixed model corres- 
ponds to just using the boundary conditions (4.2a) and letting the mixed layer 
depth equal H. Third, just as in Section 3, these equations reach equilibrium 
when surface degassing balances sparger injection. Therefore, equilibrium 
requires the same value for surface concentration as found in (3.3). However, 
with vertical diffusion tc increases with depth; thus, for a given value of tc 
at the surface, the reservoir has a higher capacity for CC^. The retention 
time increases from the estimate in Section 3, both because some degassing 
occurs as the CO 2 concentration rises and because the basin capacity 
increases. 

We illustrate the character of the time dependent solutions in 
Figure 7 which shows the rate of surface degassing relative to the sparger 
injection rate for a series of models with basin size 500 km. The three 
curves are for (1) a vertically well-mixed model, (2) K(V) = 4000 m 2 yr" 1 
which is the most reasonable estimate for the diffusion coefficient, and (3) a 
model with K(V) = 400 m 2 yr“*, an order of magnitude lower than the diffusion 
coefficient typical of the global ocean. These models all display an initial 
phase where the degassing rate increases relatively rapidly as CO 2 begins to 
fill the basin, followed later by a slow relaxation to steady state where 
degassing would equal the influx from the sparger. 

While the model takes infinite time to reach the exact steady state, 
we can assign a nominal retention time, t r , based on the condition that CO 2 
degassing equals some large fraction of the injection rate, say 90%. Figure 8 
shows the retention time estimated in that way for a series of models with 
varying radius and diffusivity. Note that vertical diffusion increases the 
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retention time relative to well mixed models, as discussed above, but for our 
best estimate of K(V), the retention time is still only about 3.5 years even 
for the largest basin size. 

Figure 9 shows the variation with depth of CO 2 concentration tc and 
pH for models in steady state. In steady state the flux of CO 2 is constant 
with depth so K(V)dtc/dz = F$ A . Thus the C0 2 concentration profiles are 
nearly linear with depth. (The relation becomes curved inside the bubble 
plume where the source function varies with depth.) Notice that models with 
longer retention time also have higher concentrations of CO 2 at the seafloor. 

Finally, we use this model to illustrate sensitivity of retention 
time to sparger depth. Figure 10 shows a series of models with K(V) = 

4000 m 2 yr"*, our best estimate, and a range of values for H, the sparger 
depth. In each case, we include a bubble plume of 30 m and a mixed layer of 
30 m, but we vary the depth H from 60 to 1000 m. The results follow the basic 
scaling appropriate to propagation of a diffusion front over a distance L 

L 2 « Kt (4.3) 

Thus, retention times become dramatically longer as the sparger depth 
increases . 

Our results showing very long retention times for disposal at great 
depth are in agreement with those of Hoffert et al . (1979) who investigated 
disposal in the deep ocean at 4000 m. If currents along the bottom of the 
South China Sea basin fortuitously carried CO 2 from the Natuna site to deep 
water, then retention times would be increased. An area of deep water 
(>1000 m) does exist a few hundred kilometers to the northeast. Although 
currents could carry CO 2 over such distances, buoyancy effects are likely to 
prevent currents from penetrating to great depth. 


17 


b. Rectangular Models With Horizontal Advection and Vertical Diffusion 

Models in Section (4a) were shown for a series of basin sizes, but 
nothing in the physical description established the basin size. In nature the 
interplay between horizontal and vertical transport establishes a scale: CO 2 
is not lost until it encounters the surface but, during the time taken to 
rise, CO 2 also migrates horizontally. For the Natuna site there is no ques- 
tion that significant horizontal motion occurs, the problem is that time and 
velocity scales are sufficiently long that basin scale circulation must be 
included to evaluate the spreading. With a typical velocity of 30 cm s" 1 
directed southeast at depth, CO 2 injected at the Natuna site reaches the coast 
of Borneo, 500 km distant, in only 20 days, which is short compared with the 
retention time (unless rapid vertical mixing occurs). If the current persists 
that far, then at Borneo the direction will deflect along the coast, and enter 
into circulation on a basin scale size, for which three dimensional time 
dependent data do not exist. 

Bearing in mind the preceding discussion, we consider two dimen- 
sional advective, diffusive flow merely to establish the general character of 
the solution, rather than to model basin scale circulation. We consider 
uniform advection with velocity v along the x direction, in a rectangular 
basin of width W, transverse to the flow. W is not determined by the model. 

We assume that CO 2 mixes rapidly along the width, and that the CO 2 source is 
distributed smoothly along a 10 km section in x of width W. The model is 
defined by the following equations and boundary conditions 


ftc = a_ K(v) |ic . Siv^tcl + q(x>2) 

atc(M) _ „, |M ate, a[v*tc(M)] f sa 

at " [ 1 az 1 H-h M ' ax " h M 


(4.4) 


(4.5a) 
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!r = ° 

at z = 0 

(4.5b) 


at x = Xj 

(4.6a) 

9tC _ q 

9x 

at x = x^ 

(4.6b) 


where x-p x 2 denote the endpoints of the rectangular model. The boundary 
conditions in (4.6a,b) apply at the "x" boundaries of the rectangle. These 
conditions imply that there will be zero diffusive flux through the 
boundaries. This is a good approximation since advection dominates the flow. 

We consider two models of the advective diffusive class, both have 
width W = 10 km, with a sparger source distributed along 10 km length in x. 

The first model is well mixed vertically (K(V) -► «), the second has K(V) = 

4000 m 2 yr _1 . Figure 11 illustrates the variation of concentration with dis- 
tance for the well mixed model. Here the total length of the model is 10 4 km, 
with the sparger source centered at x = 505 km. Because the model is well 
mixed, degassing commences directly over the sparger, and the concentration 
decays with distance downstream from the sparger. The advective timescale, 

(x 2 - xj)/v = 1 yr, controls approach to steady state in this model, so that 
concentration achieves steady state values as the flow front passes a given 
point. At this relatively high rate of flow, even within 10 4 km of the 
sparger, the net rate of surface degassing only reaches about 50 percent of 
the sparger injection rate. Further degassing of C0 2 would occur beyond x 2 , 
the downstream limit of our computational domain. However, we terminated the 
calculation at 10 4 km because the geometrical distance greatly exceeds the 
size of the South China Sea basin. 



19 


The second advective diffusive model resembles the first, but 
includes vertical diffusion with K(V) = 4000 m 2 yr' 1 rather than being well 
mixed. The basin length in this model is only 10^ km, with the sparger 
located at 95 km. For this model both the advective and diffusive times play 
a role in the approach to steady state. Here advected flow crosses the grid 
horizontally in a time (X 2 - x^)/v = 0.1 year, but the vertical diffusion 
time, H 2 /K(V) = 3 years, is much longer. Figure 12 illustrates the time 
dependent solution for the total rate of surface outgassing over the entire 
grid. Net outflow only approaches steady state after 3 years, the diffusive 
timescale. At that time, the surface outgassing rate is 35 percent of the 
injection rate. Again, as in the previous model, further degassing would 
occur beyond the downstream boundary at X 2 » Figure 13 shows contours of 
constant pH for this model after the flow achieves equilibrium. Notice that 
pH indicates large increases in acidity near the sparger. 

The essential point to be learned from the advective diffusive 
models discussed here is that the vertical diffusion time is sufficiently long 
that advection can transport C0 2 over distances large compared with the basin 
size of the South China Sea. It is unrealistic to assume that such flow would 
remain confined to a geometry resembling flow in a basin of rather narrow 
rectangular cross section as assumed here, i.e. with W constant. It follows 
that current behavior far from the Natuna site might control the ultimate fate 
of CO 2 . For instance, it is possible that currents flowing at depth near the 
site, might upwell as they approach the coast of Borneo, so that outgassing 
occurs at that point. On the other hand, currents could carry CO 2 to deep 
water increasing the retention time. 

One other point learned from this model is that, when advection 
dominates the flow, the local concentration in the vicinity of the sparger is 
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governed by the considerations described in Section 2. There we analyzed for 
peak local concentration simply by constructing the ratio of the CO 2 injection 
rate to the water flow rate. The important point from that analysis was that 
C0 2 concentrations become very large, and pH is strongly reduced, during 
periods when the flow velocity stagnates over the sparger. 

c. Axisymmetric Models With Radial and Vertical Diffusion 

Guided by insight from the previous models we consider one final 
model in which diffusion occurs in both the vertical and radial dimensions, 
but without advection. The motivation for this model is an attempt to 
approximate the effect of variable horizontal currents on transport as a 
diffusion process. Our estimate for the horizontal diffusion coefficient K(H) 
obtained in this fashion, agrees approximately with oceanographic estimates 
for K(H) suitable for open ocean conditions. Such estimates show that the 
ratio K(H)/K(V) is of order 10 7 or more (Kuo and Veronis, 1973; Sarmiento, 
1982). The enormous enhancement occurs because horizontal dynamical processes 
are not damped by stratification, as in vertical diffusion. 

We estimate the horizontal eddy diffusivity for random horizontal 
current patterns in which the velocity persists at a given value for a 
specific length of time At, and then changes its direction completely. For 
such a model, the Fickian diffusion coefficient is the average value of the 
product of length and velocity <Lv>. To estimate a characteristic value for 
K(H) we use v = 30 cm s” 1 , and a persistence time of one day 

K(H) = <Lv> - v 2 At 

K(H) - 2.5 x 10 11 m 2 yr" 1 ( — ^ (— ) 

30 cm s” day 


(4.7) 
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We stress that we do not know the actual values of K(H) or K(V) at the Natuna 
site. Our estimate for K(H) is solely for motivational purposes. It does, 
however, agree roughly with estimates characterizing the open ocean. 

In the model below we use K(H)/K( V) = 5 x 10 7 . For this ratio, and 
a typical basin size, diffusive transport over the vertical distance, 100 m, 
occurs in times comparable with transport horizontally over 500 km. This 
follows from (4.3) since (R/H)^[K(V)/K(H)] = 0.5 

With this picture for diffusion we can model the basin as an axi- 
symetric two dimensional diffusion problem: 


If =q(r,z) a^KIH) r^tc ♦ fe K(V) ff 

(4.8) 

atc(M) 1 9 ,/un 9 K(V) 9tC, F SA 

at r 9r 9r tc h M 9z H-h M h M 

(4.9a) 

—1 „ = o 

9Z 1 Z = 0 

(4.9b) 

—1 „ = 0 

9r 1 r = 0 

(4.10a) 

—1 n = 0 

9r 1 r = R 

(4.10b) 

that the time dependent variation of outgassing 

for a single 


model with a radius of 500 km, typical of the basin size. The model reaches 
90% outgassing in 6.5 years. 

Without detailed data and simulations of three dimensional time 


dependent C0 2 transport, the model described here represents our most real- 
istic simulation of the C0 2 retention time. From previous models it is clear 
that, for typical values of vertical diffusivity, C0 2 will mix to the surface 
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in a timescale of a few years, after which outgassing to the atmosphere will 
match CO2 injection rates. Also, the models, and simple physical estimates, 
show that a timescale of a few years allows transport to distribute CO2 
throughout the South China Sea basin. 


5 . SUMMARY 

Because of concern about the CC>2 greenhouse issue, seawater sparging 
of waste gas CO2 from the Natuna project was proposed as an alternative to 
direct atmospheric release. However, our results indicate tht sparged CO2 
reaches the surface and degasses to the atmosphere within a period less than a 
decade, which is short compared with the 100 year period projected for CO2 
doubling. Furthermore, increasing the concentration of CO2 in seawater 
increases acidity. 

For our best estimates of transport properties, vertical diffusion 
from 140 m depth to surface only requires 6.5 years before the rate of CO2 
outgassing through the sea surface nearly equals the sparger injection rate. 
Those estimates are based on parameter values determined from studies of the 
global ocean, conditions could be different in the shallow basin of the South 
China Sea. If more rapid vertical mixing occurs, the outgassing occurs even 
faster. Possibly, current patterns might transport CO2 below the surface for 
rather large distances, and then upwell in localized regions, perhaps near 
coastlines. Once CO2 contacts the surface at high concentration rapid 
degassing would occur. Alternatively, it is possible that currents could 
carry CO2 to deeper water giving a longer retention time, but buoyancy effects 
probably prevent currents from flowing to great depth. 
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For typical diffusion rates, the only way to keep CO 2 from degassing 
quickly is to sparge at far greater depths. Below 1000 km the CO 2 retention 
time exceeds 100 years. A region of much deeper water is available several 
hundred km to the northeast, but sparging there would be far more costly since 
a pipeline would be required, and waste gas must be pumped to greater depth 
and pressure. 

CO 2 in seawater is not inert; increasing the concentration of CO 2 
decreases pH (increases acidity). Present ambient conditions have pH = 8.3 in 
our models. Reduction of 0.7 units to 7.6 pH allows calcium carbonate in the 
aragonite form to begin to dissolve. Our models indicate that pH reductions 
of one unit or more will affect regions of order 1000 km . A tradeoff 
occurs: if the vertical diffusion coefficient decreases, so that CO 2 is 
retained longer, then the concentration of CO 2 at depth rises. Consequently, 
when the retention time increases, the impact associated with pH change also 
increases. Either the size of the affected region grows or the magnitude of 
the pH change increases, or both. 

We also estimated the peak local concentration of CO 2 and associated 
change in pH, based on advection dominated flow in the vicinity of the 
sparger, where concentration changes maximize. Those estimates show appreci- 
able pH changes whenever flow speed falls below 50 cm s"^, which occurs 
commonly. Furthermore, changes rise dramatically during stagnation events 
when flow speed drops below 10 cm s“^ which do occur occasionally. 

Our conclusion is that sparging of C0 2 offers little advantage over 
direct atmospheric release of C 02 > since the retention times is short in any 
case, but that sparging of large amounts of CO 2 does cause a negative impact 
in seawater by its affect on pH. 
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H 

Basin depth 

130 m 

h B 

Height of bubble plume 

30 m 

h M 

Depth of surface mixed layer 

30 m 

f sa 

C0 2 flux from sea to air 

moles/(m^ yr) 


F sa = E[P s (C0 2 )-P a (C 0 2 )]/340 yatm 


E 

Exchange coefficient for F$ A 

20 moles/(m^ yr) 

K ( V ) 

Vertical diffusion coefficient 

4000 m^ yr"-*- 

K(H) 

Horizontal diffusion coefficient 

2 x 10^ m^ yr"^ 

V 

Horizontal current speed 

0.30 m s _1 
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FIGURE 1. 

FIGURE 2. 

FIGURE 3. 

FIGURE 4. 

FIGURE 5. 

FIGURE 6. 


FIGURE CAPTIONS 


Geography of the Natuna site. 

Schematic of the sparger configuration. 

Variation of P(C02) and pH as a function of alkalinity TA and 
total carbon tc (per unit volume) at a salinity of 35 o/oo and 
temperature 20°C corresponding to bottom water at the Natuna 
site. Calcium carbonate is unsaturated and dissolves when pH 
< 7.4. 

As in Figure 3, but for 30°C corresponding to surface water at the 
Natuna site. 

Variation of local total carbon abundance tc (a), P(C0 2 ) (b), and 
pH (c) with current, for flushing at a uniform flow rate over the 
10 km sparger diameter. C0 2 is injected at a rate of 1.5 x 10^ 
moles per year uniformly over the bubble column of height 30 m. 

Variation of carbon abundance TC (a), P(C0 2 ) (b), retention time i 
(c), and pH (d) as a function of basin diameter R such that a well 
mixed reservoir achieves steady state C0 2 concentration. For each 
model C0 2 degasses to the atmosphere at the same rate as the 


sparger injects C0 2 . 
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FIGURE 7. CO 2 surface degassing rate versus time for a horizontally well 

mixed cylindrical model with radius 500 km. The degassing rate is 
scaled relative to the CO 2 injection rate from the subsurface 
sparger. The top curve is for a vertically well mixed model. The 
middle curve is for our best estimate of the vertical diffusivity. 
The bottom curve is for diffusivity a factor of ten below our best 
estimate. 

FIGURE 8. Retention time as a function of basin size and diffusivity for the 
cylindrical model with vertical diffusion and rapid horizontal 
mixing. The retention time is defined as the time when surface 
degassing of CO 2 reaches 90% of the sparger injection rate. 

FIGURE 9. Variation of (a) CO 2 concentration and (b) pH with height above 
seafloor for models with varying basin size and diffusivity. 

These models have reached steady state. With diffusion the 
concentration varies linearly with depth except in the mixed 
layer, where concentration is constant, and inside the bubble 
plume layer where the source function varies. 

FIGURE 10. CO 2 surface degassing rate versus time as a function of variable 
sparger depth, H. As in Figure 7 the degassing rate is scaled 
relative to the sparger injection rate, and the model is cylin- 
drical with rapid horizontal mixing and vertical diffusion. For 
all curves K(V) = 4000 m 2 yr"*. The results show that retention 
time increases roughly as H 2 for fixed diffusivity. 
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FIGURE 11. Time dependent solution for a well mixed model with advective 
velocity 30 cm s“* and width 10 km: (a) carbon concentration, 

(b) surface flux F<^, (c) total surface CO 2 outgassing relative to 
sparger injection. The 10 km long sparger source is centered at 
505 km. 


FIGURE 12. Variation of the surface degassing rate with time for a two dimen- 
sional advective, diffusive model. The model is 1000 km long by 
10 km wide and 140 m high, with the sparger source centered at 
95 km. The horizontal advection velocity is 30 cm s"*, and the 
vertical diffusivity K(V) = 4000 m 2 yr"*. After 3 years the 
degassing rate within 900 km of the sparger equals approximately 
35 percent of the injection rate. 

FIGURE 13. Contours of pH for the advective diffusive model shown in 
Figure 12, after attaining steady state. 

FIGURE 14. Variation of the surface degassing rate with time for an axisym- 
metric model with both vertical and horizontal diffusion. In this 
model K(V) = 4000 m 2 yr _1 , K(H)/K(V) = 5 x 10 7 , H = 140 m, and the 


basin radius is 500 km. 
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